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ABSTRACT 

Calculated and neasured valies of helicopter rotor flapping acgles 
in forward flight are conpared for a model rotor in a wind tunnel and 
an autogiro in gliding fli^t. The lateral flapping angles can be 
accurately predicted when a calculation of the nonunifora wake> induced 
velocity is used. At low aidvauice ratios, it is also necessairy to vise 
a free wake geometry calculation. For the cases considered, the tip 
vortices in the rotor wadce remain very close to the tip>path plauie, so 
the calculated values of the flapping notion are sensitive to the fine 
details of the wake structure, specifically the viscous core radius of Uie 
tip vortices. 


*Heau) of Hotorcraft Research Section, Large Scale Aerodynaunlcs Branch 



INT90IUCTI0N 


The blades of an articulated helicopter rotor have a hinge at the 
blade root that allows out-of-plane notion called flapping. In forward 
fU^t a once-per-revolutlon notion of the blade about the flap hinge Is 
produced, that corresponds to longitudinal and lateral tilt of the rotor 
tip-path plane relative to the shaft. Such notion occurs with hingeless 
and bearingless helicopter rotor designs as well, due to blade bending at 
the root rather than rotation about a flap hinge. The slnplest amalyses 
of helicopter rotor behavior assune that the Induced velocl^ distribution 
Is constant over the rotor disk. In fact however, the Induced velocl'^ 
distribution Is highly nonunlfom, and Uie rotor blade 'lapping notion 
Is quite sensitive to the Inflow distribution. Consider a longitudinal 
gradient of the Induced velocity. Increasing fron the leauling edge to the 
trailing edge of the rotor disk. The blade lift produced by such a 
distribution will be greater on the front of the disk than on the rear, 
so the rotor will be subjected to a longitudinal aerodynamic nonent. 

An articulated rotor responds to this nonent like a gyro, so the tip-path 
plane tilts laterally, toward Uie advancing side. Similarly a lateral 
Inflow variation will produce a roll nonent on the rotor, and hence a 
longitudinal tip-path plane tilt. Ulth a hlngeless rotor, the magnitude 
of the flapping response to nonunlfom Inflow is not greatly affected, 
but the phase of the response can be significantly decreased, so that a 
longitudinal Inflow variation can produce longitudinal flapping as well 
as lateral flapping. Because of this sensitivity to the induced velocity 
distribution. It is difficult to calculate rotor blade flapping in forward 
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flight. An accurate calculation of flapping is needed in order to predict 
longitudinal and lateral cyclic control positions, particularly for new 
helicopter designs. Koreover, since the hlade flap notion represents 
the integrated effect of the l/rev aerodynanic environnent of the rotor, 
difficulties in calculating it have serious implications regaunling the 
predictive capability of rotor aerodynamic analyses in general. 

In 193^. Vheatley^ extended the theory of Glauert and Lock for the 
calculation of helicopter rotor forces and blade notion in forward flight. 

He evaluated the accuracy of the theory by conparing with fli^t test 
results for a Pitcairn autogiro in a glide, the rotor angle-of-attadc. 
rotational speed, thrust coefficient, and longitudinal flapping angle 
were predicted with good accuracy, up to an advance ratio of about 0.4 or 
0.5. The prediction of the lateral flapping angle was consistently low 
however. Wheatley identified the use of uniform inflow as the most critical 
assumption of the analysis. The sensitivi^ of the lateral flapping calculation 
to nonunifom inflow was demonstrated by using an induced velocity distribution 
that varied linearly from the leading edge to the trailing edge of the 
rotor disk. With such a longitudinal inflow variation, the prediction of 
the lateral flapping was improved, but was still significantly low (for 
an aunplitude of the variation equal to 50^ of the mean value). Wheatley 
concluded that Ihe rotor blade flap motion can not be rigorously calculated 
without an accurats determination of ttw wake-induced inflow velocities. 

2 

In 1972 . Harris conducted a wind tunnel test of a model helicopter 
rotor, to examine the influence of advance ratio, thrust, and shaft 
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angle>of-attack on the rotor flapping angles* He found a large discrepancy 
between the measured lateral flapping and the values calculated using 
uniform Inflow. Longitudinal flapping was estimated fairly well using 
uniform Inflow* as were the trends with thrust and shaft angle. By using 
a nonuniform Induced velocity calculation based on undlstorted wake geometry* 
the prediction of lateral flapping was Improved* but the measured results 
were still significantly underestimated at low advance ratio. Harris 
suggested that the primary source of the remaining discrepancy was the 
use of undlstorted wake geometry. 

The conclusions of Wheatley and Harris have remained conjectural. 

It is the purpose of this paper to examine the Influence of nonuniform 
Inflow and wake geometry on the calculation of helicopter rotor lateral 
flapping auigles. The intent also is to establish the adequacy of the 
wake model In a helicopter rotor analysis, to the exter that may be 
Inferred by comparisons with measurements of blade flapping motion. 

ANALYTICAL MOEEL 

The TOtor flapping motion was calculated using the helicopter analysis 
described in reference 3 * The wake analysis was based on a discrete element 
representation of the vortlcity* with models for the wake rollup and 
distorted wake geometry. The free wake geometry was calculated using the 
analysis of reference 4, The blade aerodynamic loading was calciJ^ated 
using lifting line theory, with corrections for three-dimensional effects 
at the tip and at blade-vortex interactions. The rotor blade flap motion 
was calculated by means of a harmonic analysis method. Rigid body flap 


- 4 - 



motion was the only degree of fMedon oonnldered* Four harmonics of the 
blade motion were calculated, using an azimuthal Increment of 15° • 

Including blade lag notion, bending motion, or up to 10 harmonics In the 
calculation had little effect on the tip-path plane tilt. The effects of 
the blade torsion motion are discussed later In the paper. The flap 
motion relative to the shsift may be written as a Fourier series t 

^ + ^j^^cosHi + ^j^gSln*p + higher harmonics 

Here Is the longitudinal tilt angle of the tip-path plane relative 
to the shaft, positive for forward tilt; Is the lateral angle, positive 

for tilt toweurd the retreating side of the disk; and Is the blade 
azimuth angle, measured from downstream. If there Is no cyclic pitch, 
both and will be negative In forward flight. 

The Induced velocity, circulation, and lift were evaluated at 
15 stations along the blade, concentrated toward the tip. The trailed 
■ vortlclty In the wake directly behlxul the blade was represented by discrete 
vortex lines positioned midway between the points at which the circulation 
was calculated. The strength of each trailed line was defined by the 
difference between the bound circulation at successive radlzd stations. 

This part of the model Is a common numerical Implementation of lifting 
line theory. For the rotary wing It Is necessary also to model the rollup 
of the vortlclty Into a concentrated tip vortex, because the blade 
encounters the wake from the preceding blades as It rotates. Therefore, 
after an azimuthal extent of 30° In this case, the trailed vortlclty was 
concentrated into a single tip vortex line, with strength equal to the 
maximum bound circulation of the blade. For conservation of vortlclty. 
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there must be an Inboard sheet of trailed vortleity with equal total 
strength and opposite sign as the tip vortex* This Inboard sheet Is much 
less Important that the tip vortex, partly because the vortleity Is not as 
concentrated and partly because It tends to be convected downward faster 
than the tip vortex. Consequently the entire Inboard sheet could be 
reasonably represented by a single vortex line, with a lairge core radius 
(about of the rotor radius In this case) to avoid unrealistically large 
Induced velocities when the wake passed under the following blade. The 
bound circulation In forward flight varies azlmuthally as well as radially, 
so the shed wake must also be Included. The shed wake was modelled by 
radial, discrete vortex lines (also with a large core radius). A similar 
wadee model was used for the other blades of the rotor. The detailed near 
wake model need not be used for these other blades since It Is fair from 
where the Induced velocity is being calculated on the first blade. The 
curved trailed vortex elements in the wake are represented by a connected 
series of straight line segments, with an azimuthal increment of 15° In 
this case. The vortleity strength varied llneaurly along each vortex line 
element, both trailed and shed. Four revolutions of the wake behind each 
blade were modelled in this fashion. It was determined that three revolutions 
were sufficient, except at the very lowest advance ratio analyzed. 

A distributed vortleity model was used for the tip vortex core, 
in which one-half the vortleity was outside the core radius (defined at 
the point of maximum tangential velocity). With this model, the maximum 
velocity Induced by the vortex was half that of a vortex with all the vortleity 
concentrated inside the core radius. 
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An undlstorted wake geometry Is obtained by assuming that each 
element of vorticlty in the wake is convected downward at a rate equal to 
the mean induced velocity over the rotor disk, as the blades rotate and the 
helicopter flies forward. The resudting geometry consists of skewed, 
interlocking helices, one behind each blade. The self>induced distortion 
of the wake is neglected in this model. The analysis of reference 4 weis used 
to calculate the distorted, free wadce geometry in forward flight. The 
free wake geometry was calculated for the tip vortices only, for two 
revolutions of the wake behind each blade. Only two iterations within 
the free wake geometry calculation were necessary. 

The analysis progressed in three stages. In the first stage, the 
blade motion and loading were calculated using uniform inflow. Next the 
wake influence coefficients were calculated with an undistorted wake geometry, 
and the motion and loading obtained using nonuniform inflow. Finally, 
the free wake geometry was calculated, the influence coefficients were 
re-evaluated, and the blade motion and loading again obtained using nonuniform 
inflow. It was determined that omitting any of these three stages gave 
unconverged results, but it was not necessaury to repeat emy stage before 
going on to the next. Within each stage, ihere was an iteration so that 
the calculated induced velocity and bound circulation were consistent. 

Also wiUiin each stage, there was an Iteration on the collective pitch, 
in order to trim the rotor thrust to a prescribed value. 

COMPARISON WITH WIND TUNNEL DATA 

Harris measured the flapping angles of a model rotor in a wind tunnel. 

The principal parameters defining the rotor and operating conditions are 
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Table 1. Rotcr Ibrameters and Operating Conditions 



Model Rotor 

Autogiro Rotor 

Radius, R 

0.832 m 

6.86 m 

Number of blades 

k 

4 

Solidity, e- 

0.0891 

0.1037 

Lock number 

5. 80 

19.1 

Flap hinge offset, e/R 

0.0229 

0.0125 

TWist 

-9.1U® 

5.° 

Longitudinal cyclic 

0.73° 

0° 

Lateral cyclic 

0° 

0° 

Tip speed, CLr 

137. m/sec 

102. m/sec 

Nominal C,j/«- 

0.080 

0.064 

Airfoil 

V23OIO-I.58 

Gottingen 429 

Source of data 

Reference 2 

Reference 1 
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given in table 1. Complete details of the rotor and test data may be found 
in reference Z, An accuracy of db0.25° is given for the flapping measurements* 
The wind tunnel had partly slotted walls* designed to produce no angle»of> 
attack change due to wall constraints* and the rotor diameter wc ' o.:ly 
Z7fo of the tunnel width. Flow breakdown was estimated to occur at advance 
ratios below 0.04, 

Figure 1 compares the calculated and measured values of the lateral 
flapping angle as a function of advance ratio* for a thrust coefficient 
to solidity ratio of 0.08 and a tip* path plane angle*of*attadc of approximately 
1° (aft tilt). For reference* a helicopter at an advance ratio of 0.13 
would typically have a tip-path plane angle-of*attadi of approximately 
forward. The lateral flapping is underpredicted when uinlform inflow is 
used, euid even when nonuniform inflow based on the undistorted wake geometry 
is used. Below an advance ratio of about 0.l6* it is necessary to include 
the free wake calculation in order to obtain a good estimate of the lateral 
flapping. Figure Z shows the corresponding results for the longitudinal 
flapping angle, for which there is much less influence of the inflow model. 

The free wake geometry analysis was not applicable at hover and very low 
speed, so no calculated results are shown for advance rat ..os below 0.04. 

It was also found that the induced power calculation was significantly 
influenced by the wake model. Below an advance ratio of about 0.10, the 
calculated induced power was actually less than the ideal momentum theory 
value when the xindistorted wake geometry was used* but was appropriately 
greater than the ideal value with the free wake geometry. For advance 
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ratios greater than 0.15i the calculated induced power was about the same 
for the two models. 

Figure 3 shows the calculated vertical displacement of the tip 
vortices relative to the tip-path plane, for the two weJce geometry models. 

The wake geometry is shown as a function of the wake age <j> (the azimuth 
angle along the wake helix), for the rotor blade at four azimuth angles. 

The significant self- induced distortion of the free wake ge. metry is evident, 
resulting in numerous blade-vortex interactions in which the vertical 
separation is a fraction of the blade chord (which is 7^ of the rotor radius). 
The tip vortices even pass above the blades initially, a behavior that has 
been observed experimentally as well.^ Figure 4 shows the influence of 
the wake geometry on the calculated induced velocity. The discrete tip 
/ortices of the rotor wake produce a l/rev variation of the induced 
velvcJty that corresponds to a longitudinal gradient over the rotor disk, 
and hence increases the lateral flapping. By moving the wake closer 
the tip-path plane, the induced velocity change on the sides of the d:. 
where the blaides sweep over the tip vortices, is greatly increased. 
Consequently the wake geometry is an important facto, in detnnnining the 
lateral flapping at low advance ratio and low tip-path plane tilt. 

In the cases considered here, the free wake geometry places the 
tip vortices so close to the blades that the calculated flapping is sensitive 
to the value of the viscous core radius, which determines the maximum 
velocity induced by the vortex. For the calculations presented in figures 
1 to 4, a core radius of 0.05R was used. No measurements of the core 
radius are available for the model rotor being analyzed, so this value 
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was chosen based on the correlation produced. Figure 3 shows the sensitivity 
of the calculated Induced velocity to the tip vortex core radius, and 
table 2 shows the corresponding calculated flapping angles. With the 
undistorted wake geometry, there is little Influence of the core radius. 

These values of the vortex core radius should not be Interpreted too 
literally however. There are a number of factors in addition to the 
core radius that combine to determine the magnitude of the vortex- Induced 
loading, including the tip vortex stength, the extent of the tip vortex 
rollup, lifting surface effects on the induced blade loading, and possibly 
even vortex bursting or vortex- induced stall on the blade. In the absense 
of complete information about each of these phenomena, the vortex core 
radius is a convenient parameter with which to account for their cumulative 
influence on the rotor blade loading. 

Figure 6 compares the measured and calculated values of the lateral 
flapping as a function of thrust, at an advance ratio of 0.08 and shaift 
angle of -1.33*^« 'The free wake geometry is required to accurately predict 
the rate of change of with C,jyV-, and it is Increasingly important at 
high thrust. Figure 7 compares the measured and calculated values of the 
lateral flapping as a function of the shaft angle-of-attack, at an advance 
ratio of 0.08 and a collective pitch of 11°. 

In all cases, the collective pitch was adjusted to trim the calculated 
thrust to the measured value. It was found that the calculated collective 
pitch angles were consistently lower than the measured values, a common 
result in rotor analyses that may usually be attributed to elastic pitch 
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T^ble 2* Calculated Influence of vortex core radius on blade 
flapping angles “ 0.1, - 0.08, and 

Flapping Angles (deg) 

Pic Pis 

Free Hake geometry 



core radius, r O.OIOR 

c 

-2.08 

-b.l7 

0.030R 

-2.b7 

-3.39 

0.050R 

- 2.57 

- 3.13 

0.075R 

- 2.56 

-2.t»7 

O.IOOR 

- 2.52 

-2.68 

Undistorted wake geometry 



core radius, r » O.OIOR 
c 

-2.77 

-1.69 

0.Q30R 

-2.73 

- 1.70 

O.O 5 OR 

-2.70 

- 1.70 

0.075R 

-2.66 

-1.68 

O.IOOR 

- 2.62 

- 1.65 
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deflection of the blade* Insufficient information was available regarding 
the physical characteristics determining the blade torsion dynamics (control 
system auad blade stiffness* blade inertia and eenter-of>gravity distributions* 
and the aerodynamic center distribution) to alloN an exi^oration of this 
facet of the analysis. 


OOHPARISOM WITH PLIGHT SdTA 

Wheatley measured the flapping angles of an autoglro in gliding 
flight. The rotor and operating parameters are summarised in table it 
complete details of the aircraft and test data are given in references 1 
and 6. An accuracy of ^ 0*1^ is given for the flapping measurements. 

In the calculations the rotor was trimmed by setting the shaft angle-of- 
attack to the measured value* and adjusting the collective pitch to obtain 
the measured value of the rotor thrust coefficient. On an autogiro the 

rotor has a positive aogle-of-attadc (aft tilt of the tlp>path plane)* 
that becomes fairly large at low advance ratio. 

Figure 8 compares the calculated and measured values of the lateral 
flapping as a function of advance ratio. Again it is essential to use 
nonuniform inflow in order to accurately predict the flapping angles* 
and below an advance ratio of about 0.20 it is also necessary to use the 
free wake geometry. Figure 9 shows the corresponding results for the 
longitudinal flapping* for which the influence of the nonunifom inflow is 
smaller but not negligible* 
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OOMCXUUDiC REMARKS 


Cadculated and neasumd values of helicopter rotor flapping angles 
in forward flight have been conpeured for a aodel rotor in a wind tunnel 
and an autogiro rotor in gliding flight. It has been confined that the 
discrepancy between the experlsental data and theories based on uniform 
inflow is attributable prinarlly to the highly nonunifon velocity indticed 
at the blades by the tip vortices in the rotor wake. By using an analysis 
based on a calculation of the nonunifon inflow* it was possible to accuntely 
CEilculate the rotor lateral flap^dog angles. At low advance ratios it is 
also necessary to use a free wake geoeetry. a result anticipated hy Harris 
although not by Wheatley, the degree of sensitivity of an integrated effect 
of the rotor aerodynaaic environnent (the lateocal flapping an^e) to the 
fine structure of the wake is perhaps surprising, but plausible upon 
exaainatlon of the phenoaenon. It is concluded that the principal features 
of the rotor wake are correctly represented in the nonunifon inflow and 
free wake geonetry aodels utilized in this investigation. To correlate 
and develop these models further will require experimentakl data that 
encompasses the fine details of the rotor aerodynamics, such as measirements 
of the tip vortex circulation, core size, and peak velocities, even when 
the parameters of most interest are the integrated effects such as the first 
harmonic flapping motion. 
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Figure 1. Comparison of mesisured and calculated model 
rotor lateral flapping angles as a function 
of advance ratio • 0.08 and 
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O EXPERIMENT 
UNIFORM INFLOW 

— — NONUNIFORM INFLOW. UNOISTORTEO WAKE 
— — NONUNIFORM INFLOW. FREE WAKE 



Figure 2. Conparison of measured and calculated model rotor 

longitudinal flapping angles as a function of advance 
ratio (dj/r- » 0.08 and «<^pp ~ 1°) 
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f GENERATING BLADE 

- ■ - FREEWARE 

— — UNDISTORTED WAKE 
TIP-PATH PLANE 

0 BLADE-VORTEX INTERACTION 



. Calculated tip vortex vertical dlsplacemeat 
(yU- 0,1, C,j/o- “ 0.08, and 1°) 
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Flgxire 3 



UNIFORM INFLOW 

NONUNIFORM INFLOW. UNOISTORTED WAKE 

NONUNIFORM INFLOW. FREE WAKE 



Figure 4. Calculated induced velocity at ladiaj. station 
r/R “ 0.92 for the three wake aiodels ( JK » 0.1, 
Oj/e- - 0.08, and « 1°) 
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rc/R 

0.01 



vi'.deg 

Figure 3 * Calculated induced velocity at radial station 

r/R = 0«92 tasing the free wake model, as a function 
of tip vortex core radius 0«1* C,j/v - 0.08, 

and 1°) 
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O EXPERIMENT 


UNIFORM INFLOW 

NONUNIFORM INFLOW, UNDISTORTED WAKE 
NONUNIFORM INFLOW, FREE WAKE 



rotor lateral flapping angles as a function 
of thrust “ 0.08 and = - 1.35°) 


- 21 - 


o 


EXPERIMHNT 
UNIFORM INFLOW 

NONUNIFORM INFLOW, UNDISTORTED WAKE 
NONUNIFORM INFLOW, FREE WAKE 



Figure 7. Comparison of measured and calculated 
model rotor lateral flapping angles as 
a function of shaft angle » 0,08 
and ©75 - 11°) 
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O EXPERIMENT 
UNIFORM INFLOW 

NONUNIFORM INFLOW. UNOISTORTED WAKE 

NONUNIFORM INFLOW. FREE WAKE 



Figure 8. Comparison of me .sured and calculated autoglro 
rotor lateral flapping angles as a function of 
advance ratio 
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